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Abstract 


A  detailed  analysis  is  performed  for  a  finite  element  method  applied  to 
the  general  one-dimensional  convection  diffusion  problem.  Piecewise  poly¬ 
nomials  are  used  for  the  trial  space.  The  test  space  is  formed  by  locally 
projecting  L-spline  basis  functions  onto  "upwinded"  polynomials.  The  error 
is  measured  in  the  mesh  dependent  norm.  The  method  is  proven  to  be 
quasi-optimal  (yielding  nearly  the  best  approximation  from  the  trial  space) , 
provided  that  the  input  data  is  piecewise  smooth.  This  assumption  is  usually 
observed  in  practice.  These  results  are  used  to  establish  a  posteriori  error 
estimates  and  an  adaptive  mesh  refinement  strategy  in  Part  II  of  this  series 
(35). 
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CHAPTER  1 
INTRODUCTION 

Reliable  numerical  solutions  to  singularly  perturbed  boundary  value  problems 
are  of  great  importance  in  engineering.  Because  of  the  degenerative  nature  of 
these  problems,  conventional  numerical  methods  produce  approximations  whose 
optimality  degenerates  as  well. 

In  this  paper  we  consider  the  model  problem 
(1.1)  -Eu"  +  a(x)u'  +  b(x)u  =  f(x),  in  (0,1), 

u(0)  =  a, 

Sju’d)  +  B2u(l)  =  B, 

where 

u  =  u(x)  is  the  solution  and  may  measure,  for  example,  temperature  or  concentration, 
e  >  0  is  the  diffusivity  of  u(x),  a(x)  _>  a  >  0  is  the  velocity  of  the  medium 
carrying  u(x) ,  b(x)  is  the  coefficient  for  u(x)  used  to  represent  a  "loss"  if 
b(x)  >  0,  or  a  "source"  if  b(x)  <  0,  and  f(x)  is  the  external  source  term.  The 
order  of  this  problem  degenerates  from  two  to  one  as  E  +  0,  and  particular 
interest  will  be  focused  on  the  case  when  e  is  small. 

If  a(x)  !i  >  0  and  the  ratio  _a/e  is  large-in  practice  this  ratio  may  be  as 

q 

large  as  10  ,  the  solution  to  (1.1)  will  often  exhibit  boundary  layer  behavior. 

A  boundary  layer  can  be  loosely  defined  as  a  small  region  near  the  boundary  where 
the  solution  changes  rapidly.  If  the  function  f  is  "rough",  the  solution  may 
also  have  interior  layers  where  the  solution  changes  rapidly  near  some  points 
{x*}  C  (0,1).  In  cases  when  a,  b  and  f  are  smooth,  the  solution  will  generally 
behave  like  ea^^x  near  x  =  1  (see  e.g.  [  1 0] ,  [15]  or  [36]). 
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Much  attention  has  recently  been  focused  on  the  application  of  finite  element 
methods  to  problems  of  chis  type.  Because  conventional  finite  element  methods 
employ  piecewise  polynomial  trial  and  test  spaces,  one  problem  is  the  ability  of 
these  polynomials  to  approximate  the  exact  solution  well,  particularly  in  the 
boundary  layer  region.  However,  the  most  serious  problem  by  far  is  the  loss  of 
stability  or  quasi-optimality  with  these  conventional  methods. 

A  finite  element  solution  u^  is  called  quasi-opt imal  in  the  norm  11*11  if  there 
exists  a  C,  independent  of  e  and  the  mesh,  such  that 

(1.2)  | |u-uhj  |  <  C  inf  | |u-w| | , 

wOh 

where  is  the  trial  space  from  which  the  approximation  u^  to  the  exact  solution 
u  is  taken.  Whenever  (1.2)  holds  we  are  guaranteed  to  have  nearly  the  best 
approximation  from  the  trial  space  S^.  With  conventional  methods,  however,  the 
constant  C  in  (1.2)  becomes  unbounded  as  £  -*•  0. 

This  loss  of  quasi-optimality  can  be  seen  numerically  by  considering  the 
following  problem: 


(1.3a) 


-eu"  +  u’  =  f  in  (0,1) 
u(0)  =  u(l)  =  0, 


with 

(1.3b) 


f(x) 


1  if  0  <  x  <  1/3, 

0  if  1/3  <  x  <  1,  x  ^  2/3, 

6(x-2/3)  if  x  =  2/3, 


where  6(x-2/3)  is  the  Dirac  delta  function  representing  a  point  source  at  x  =  2/3. 


Besides  having  a  boundary  layer  at  x  =  1,  this  problem  has  an  interior  layer  at 
x  =  2/3.  The  exact  and  conventional  piecewise  linear  finite  element  solutions  are 
shown  in  Figure  1.1  with  £  =  .0001,  and  N  =  24  elements.  The  loss  of  quasi¬ 
optimality  expresses  itself  in  the  form  of  spurious  oscillations  of  the  finite 
element  approximation. 
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The  most  common  way  to  alleviate  the  problem  of  oscillations  is  to  "upwind" 
the  test  space.  In  the  case  of  linear  elements  this  is  done  by  adding  a 
quadratic  term,  multiplied  by  some  parameter  ct,  to  each  linear  basis  function  of 
the  test  space  (see  [  8|,  [  9],  117]-[21],  [30],  or  ,  [37]).  This  procedure 
will  be  referred  to  as  cr-quadratic  upwinding  throughout  this  paper.  A  typical 
basis  function  upwinded  in  this  way  is  displayed  in  Figure  1.2. 

In  all  of  the  papers  mentioned  in  the  preceding  paragraph,  the  criteria  used 
for  the  selection  of  «  was  either  to  eliminate  oscillations,  or  to  produce  exact 
nodal  solutions  for  the  model  problem  (1.3a).  For  example,  Christie,  et.  al.  [8], 
Heinrich,  et.  al.  [19]-[21],  Mitchell,  et.  al.  [30]  and  Zienkiewicz,  et.  al.  [37] 
have  displayed  the  "optimal"  «  which  produces  the  exact  nodal  values  for  problem 
(1.3a),  when  f (x)  =  1. 

It  is  pointed  out  by  Gresho  and  Lee  [16]  that  "ad  hoc"  upwinding  can  be 
deceptive  to  the  analyst  by  smoothing  out  the  results,  and  any  solution  obtained 
by  upwinding  does  not  represent  a  solution  to  the  original  problem.  Instead  of 
upwinding,  they  advocate  the  use  of  conventional  finite  element  methods,  and 
propose  to  use  the  information  given  by  the  oscillations  to  refine  and/or  relocate 
the  mesh  points  in  the  areas  where  any  "wiggles"  occur. 

Although  we  agree  that  the  upwinding  criteria  of  damping  the  oscillations 
is  incorrect,  the  oscillations  themselves  may  be  misleading  in  determining  how 
the  mesh  refinement  should  proceed.  The  approximation  to  problem  (1.3a,b) 
displayed  in  Figure  1.1  would  mislead  one  into  refining  the  mesh  everywhere. 

The  criteria  for  a-quadratic  upwinding  should  not  be  to  eliminate 
oscillations  but  should  be  to  obtain  quasi-optimality  as  in  (1.2)  with  the 
constant  C  independent  of  e,  and  the  mesh  spacing  h.  In  [18],  Griffiths  and 
Lorenz  attempted  to  select  a  in  a  way  to  minimize  this  constant  C.  Unfortunately, 
even  with  a  chosen  in  this  way,  if  e  <<  h  their  value  of  C  increased  with  rate 
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h-^  as  h  +  0. 

In  [  3 ] •  we  have  proven  that  quasi-optimality  is  attainable  for  problem  (1.3a) 
(in  a  norm  different  from  the  one  studied  by  Griffiths)  with  o-quadratic 
upwinding  if  and  only  if  the  input  data  f  is  piecewise  smooth  -  a  reasonable 
assumption  in  practice. 

Upwinding  can  also  be  done  with  the  use  of  L-spline  basis  functions.  Methods 
using  these  spaces  have  been  studied  by  Heraker  and  De  Groen  ([11],  [12],  [22]) 
who  prove  a-priori  estimates  of  the  error  at  the  nodes  and  in  the  norm 


In  this  norm,  however,  the  error  arising  from  any  piecewise  linear  approximation 
cannot  be  made  small  unless  h  <  e.  They  also  propose  to  upwind  the  trial  space 
in  order  to  get  a  better  fit  to  the  exact  solution  in  the  boundary  layer  region. 
However,  this  upwinding  can  introduce  spurious  internal  layers  in  the  approximation. 

Also,  whenever  the  norm  (1.4)  is  used,  the  assumption  a(x)  >  >  0  and  the 

additional  assumptions 

(1.5)  b(x)  >_  0,  and  b(x)  -  ^  a' (x)  y  Y  >  0, 

are  needed  to  prove  coercivity  of  the  bilinear  form  used  to  pose  (1.1) 
variationally.  When  a(x)  _>  a  >  0,  a  much  weaker  additional  assumption,  namely  - 

(1.6)  if  b(x)  >  b  then  a^  +  4eb  >  Y  >  0, 

is  sufficient  to  guarantee  that  zero  is  not  an  eigenvalue  for  (1.1)  and  hence  the 
solution  will  be  unique.  This  assumption  together  with  some  smoothness  of  a,  b 
and  f,  are  sufficient  for  the  results  in  this  paper.  Note  that  in  (1.6)  b  and 
hence  b(x)  may  be  negative.  Furthermore,  no  condition  on  the  expression 
b(x)  -  y  a'(x)  is  required. 

Since  the  L-spline  basis  functions  are  exponential  and  have  boundary  layers 
themselves,  a  special  quadrature  rule  must  be  devised  in  order  to  perform  the 
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integrations  needed  to  assemble  the  matrix  equations  with  sufficient  accuracy. 
Diaz-Munio  and  Wellford  [  13]  use  exponent i all v  upwinded  basis  functions  which  are 
local  asymptotic  expansions  of  the  solution,  and  describe  a  special  numerical 
quadrature  rule  which  is  exact  for  integrands  of  the  form  tneat . 

Kellogg  and  Han  [27]  present  a  scheme  in  which  they  add  one  singular  function 
of  boundary'  layer  type  to  both  the  test  and  trial  spaces.  Using  this  method  they 
were  able  to  prove  the  error  estimate 

(i- 7)  Mu  -  uh!  !1  e  <  ch, 

where  C  is  bounded  independently  of  c  and  h,  and  !  |  '  |  j ^  r  is  defined  in  (1.4). 

This  method  presumes  a-priori  knowledge  of  the  location  of  the  boundary  layer. 

For  example,  in  order  to  solve  problem  (1.3)  which  also  has  an  interior  layer  at 
x  =  2/3,  some  modifications  must  be  made  in  their  algorithm. 

In  this  paper  we  develop  a  finite  element  method  which  produces  a  quasi-opt imal 
approximation  to  (1.1)  for  all  values  of  e  C(0,ll.  The  norm  used  to  measure 
the  error  is  closelv  related  to  the  Lp  norm.  An  Lp  (and  in  particular  an  Lj) 
type  norm  is  appropriate  especially  when  the  location  of  the  boundary  or  interior 
layers  are  of  importance.  In  the  second  paper  of  this  series:  Part  II  - 
A-Posteriori  Error  Estimates  and  Adaptivity,  the  numerical  results  presented  are 
based  on  the  type  norm. 

The  norms,  spaces,  and  bilinear  form  used  to  pose  (1.1)  variat ionally  are 
presented  in  Chanter  2.  Chapters  3  and  4  show  that  a  basis  for  an  exponentially 
upwinded  test  space  can  be  found,  which  produces  a  quasi-optimal  approximation. 

In  Chapter  5  these  exponential  basis  functions  are  projected  onto  a  polvnomially 
upwinded  test  space.  This  type  of  upwinding  is  a  generalization  of  cx-quadratic 
upwinding.  Finally,  in  Chapter  6  some  remarks  are  made  on  the  extension  of  these 
results  to  more  general  problems. 
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CHAPTER  2 

MATHEMATICAL  FRAMEWORK 

This  chapter  sets  up  the  mathematical  framework  in  which  the  convection- 
diffusion  problem  (1.1)  will  be  studied.  First,  two  theorems  are  presented  which 
are  used  to  prove  existence  for  var iat ionally  posed  problems,  and  quasi-opt irnalit 
for  finite  dimensional  approximations.  Next  some  results  cone  "ning  the  Green's 
function  to  (1.1)  are  proven.  After  the  norms  and  spaces  nee  to  pose  (1.1) 
variationally  are  provided,  the  Green's  function  results  are  1  to  prove  some 
important  embedding  theorems. 

SOME  ABSTRACT  RESULTS 

Two  crucial  results  concerning  variationally  formulated  boundary  value 
problems  and  finite  element  approximations  are  given  in  this  section. 

Theorem  2.1.  Let  ^  and  ^  be  two  reflexive  Banach  spaces,  indexed  by  a 

parameter  A  with  A  varying  over  some  index  set,  with  norms  | | • | |  ^  ^  and  | | • | ^ 
respectively,  and  let  B,  be  a  bilinear  form  on  K1  .  x  K  ..  We  suppose  the 
following  are  satisfied: 


(2.1) 


I  B^ (u, v) 


±Cll 


1,A'  |v|  I  2, A  for  a11  uC  K1  ,A’  vCK2,A’ 


(2.2) 


and 

(2.3) 


inf 

uCK 


1,A 
1,  A 


=  1 


sup 

vCK 


2, A 
2, A 


I  B^(u, v)  |  >^  C2  >  0, 


1 


sup  |B.(u,v) |  >  0,  for  each  0  ^  v  C  K 
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where  C.  and  C?  are  positive  constants,  possibly  depending  on  A .  Then  if 


f  C  (K.,  . )',  there  exists  a  unique  solution  u  C  K  .to  the  problem 
2 ,  A  1  ,  A 


Moreover,  u  satisfies 


B . (u, v)  =  f  ( v) 


ul I l,i  -  C21 h f ! I K;  , 


for  each 


v  C  K2,A  • 


If  the  bilinear  form  B,(.,.)  satisfies  the  assumptions  (2.1),  (2.2)  and  (2.3), 

B.  is  said  to  be  a  (C.,  C~)-proper  bilinear  form  over  the  space  Kn  ,  x  K  .  .  It 

A  XX  X  )  Li  2  y  Li 

should  be  noted  that  (2.2)  and  (2.3)  can  be  shown  to  be  equivalent  tc 


(2.2)* 

inf 

sup 

VCK2,i 

nCK  a 

III'  =  1 

1  1  l  l  =  1 

1  1 V '  '  2 ,  A 

1 ,u' 'i.i 

and 

(2 . 3)  * 

sup  |  B  .  (u,  v) 

I  B a ( u , v ) I  >  C,  >  0, 


0  4  u  C  K-  . 

1,  A 


This  observation  will  be  specifically  used  in  tnis  paper. 

Since  we  will  be  studying  finite  element  approximations  * o  u,  we  let  ^  and 
^2  A  dimensional  subspaces  of  ^  and  ^  respectively.  Clearly, 

condition  (2.1)  holds  on  S,  ^  x  .  with  the  same  constant  C..  .  We  will  be 
invoking  the  following  theorem  concerning  the  finite  element  solution  u^. 

Theorem  2.2.  Suppose  B^  is  (Cj,C^)  -  proper  over  S^  ^  x  S0  ^  furnished  with  norms 
||*||  .  and  j  j • ]  L  A ,  respectively.  Let  u  C  K-  .,  and  let  u  C  S  .  be  the 

unique  solution  to  B^(u^,v)  =  B^(u,v)  for  all  v  C  ^  • 

Then 


-  Uh  ' 1 , A  — 


<  (1  +  7T>  inf  | |u  -  w| 


J2  wC  S, 


1,A  • 
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For  the  proof  of  Theorems  2.1  and  2.2,  see  e.g.  (  1  ]. 

We  remark  that  the  reflexivity  of  the  space  A  is  not  necessary  for 

Theorem  2.2  to  be  valid  -  the  reflcxivitv  is  needed  onlv  on  the  spaces  S,  .  and 

1,A 

S>2  ^  ,  which  is  guaranteed  since  they  are  finite  dimensional.  However, 

Theorem  2.1  will  not  imply  the  existence  of  a  solution  u£  if  ^  and 

are  not  reflexive.  Therefore,  if  K.  .  and  K0  ,  are  not  reflexive,  the  solution 
u  must  be  assumed  or  shown  to  exist  in  ^  by  some  other  method,  in  order  to 
apply  Theorem  2.2. 

THE  GREEN'S  FUNCTION 

Consider  the  operator  L  defined  by  Lu  =  -eu"  +  a(x)u'  +  b(x)  u,  for 
u£C~[0,l],  where  C^[0,1]  denotes  the  space  of  functions  with  k  continuous 
derivatives  on  [0,1].  Let  G(xjy)  denote  the  classical  Green's  function  for  the 
operator  L  with  boundary  conditions: 

G  ( 0 1  v )  =  0 

(2.4) 

r(C(x|v))  =  BjC  (1; v)  +  B2  G(ljy)  =  0  . 


The  following  lemma  will  be  used  to  establish  the  existence  of  this  function. 

Lemma  2.3.  Assume  a(x)£C^[0,l]  and  b(x)£C° [0, 1 ] .  If  there  exists  a  positive 
function  w££^[0,l]  which  satisfies 

Lw  >  0, 

(2.5)  w(0)  >  0, 

w  o, 

then  the  Green's  function  G(x|y)  to  L  exists,  is  unique,  and  is  non-negative. 
Proof:  This  lemma  follows  from  the  maximum  principle  (see  e.  g.  [32],  Chapter  1, 
Theorem  1.1),  and  standard  results  concerning  the  Green's  function  (see  e.g.  [33], 
Sections  1.3  and  1.5).  A  detailed  proof  of  this  lemma  can  be  found  in  [34], 

(Lemma  2.5). 


9 

UNCLASSIFIED 


UNCLASSIFIED 


In  order  to  find  a  function  w  satisfying  (2.5)  we  make  the  following 
assumption  for  the  operator  L. 

Al:  a(x)C  J|0,1],  a(x)  a  >  0 

b(x)  £  C°[0, 1] ,  b(x)  _>  b  and 

]j  is  such  that  _a~  +  4eb  =  Y  >  0. 

For  the  boundary  operator  T  we  assume 

A2:  Br  B2  >  0,  Bx  +  B2  >  0. 

The  following  is  a  corollary  to  Lemma  2.3. 

Corollary  2.4.  Suppose  assumptions  Al  and  A2  hold.  Then  the  Green's  function 
to  L,  satisfying  (2.4),  exists,  is  unique,  and  is  non-negative. 

Proof:  Let  a  =  (si  +  (a^  +  4e_b)  J)  >  0  and  w(x)  =  eax.  Then  w  satisfies  (2.5). 

An  important  fact  used  in  this  paper  is  that  the  Green's  function  is  bounded 
independently  of  C.  The  proof  of  the  following  theorem  which  establishes  this 
fact  is  similar  to  one  found  in  Lorenz  [29] . 

Theorem  2.5.  Suppose  assumptions  Al  and  A2  hold.  Then  the  Green's  function 
G(x|y)  for  L,  which  satisfies  (2.4),  is  bounded  by  a  constant  which  is 
independent  of  x,  y  and  e. 

Proof:  Let  G(x|y)  be  the  Green's  function  for  L  where 

A 

Lw  =  -ew"  +  a(x)w', 

with  boundary  conditions 

w(0)  =  0 

Bjw'U)  +  B2w(1)  =  0. 

Then  G(x|y)  is  given  by 
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(2.6) 


C.(x|v)  = 


BjP(l)  J~  &2  f  P (t-)dt  1  x 


ep(y)x 


p(t)dt,  x  <  y. 


Va)/op(t>dt  .  82 /^<t)dt  r1  .  . , 
- Ep(vTx  '  +  . VV(W_7  p(t)dt- 


ep(y)X 


where  p(t)  =  exp(—  j  a(s)ds),  and  >;  =  B^p(l)  +  $2  I  P(t)dt.  By  assumption  A 2, 

^  o 

61?  B2  2l  0,  and  B-^  +  &2  >  0,  and  hence  x  >  0.  Therefore,  from  (2.6)  we  have 


G(x|y) 


fl  P(t)dt  1  fy  if 

-  G(yly)  -  ep(y)"  =eJoexp(eJ, 


a(s)ds)dt 


,  1/  e(t-y)a/edt  <  (1  _  e-«/c)/a  ,  1/a- 

—  £•>  _  —  —  — 


If  b  =  0  then  b  (x)  >  0.  Let  w  (x)  =  G(x]y)  -  G(x|v) 


Lw^(x)  =  -b(x)G(xjy)  for  each  x  G  (0,1), 

Hence,  w  (x)  =  -  G(x|£)b(£)  G(£|y)d£.  Since  both  G  and  G  are  non-negative 

y  J  o 

by  Corollary  2.4,  w^(x)  0.  Thus,  G(xjy)  G(x|y)  £  1/a. 

If  b  <  0,  let  a  =  2^- (a  -  (_a^  +  4eb) ^^)  >  0,  and  L^w  =  e  '^(e^w)  =  -Ew"  + 

(a(x)  -  2eo)w'  +  (-£ 0^  +  a(x)o  +  b(x))w.  Since  -eo^  +  a(x)o  +  b(x)  _>  -ea^  + 

a  a  +  _b  =  0,  and  a(x)  -  2eo  _>  a  -  2 cc  =  (af  +  4cb)^^  =  >  0,  we  are  in  the 

case  when  b  =  0  for  the  operator  La.  Let  Ga(x|y)  be  the  Green's  function  for  La 

satisfying  the  boundary  conditions  G0(O|v)  =  0, 


do- 

FoGa=-  BlH#  (1I*>  +  «2  +  oB1)  Go(1^  "  °- 
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Since  0  0,  the  boundary  operator  Tn  satisfies  the  analogous  assumption  as 

1  /o 

A2  was  to  the  boundary  operator  T,  and  therefore,  G0(x|y)  <  1/Y 

G  (xjy)  is  related  to  G(x  y)  through  the  identity  O(xjv)  =  e°^X  ^O  ^x  |  y) . 

n  i/o  a  ?  1 ,  ~  i 

Therefore,  G(x|y)  <_  e  /y  .  Now,  0  =  — (1  -  (1  +  4eb/a~)  ")  irr  (1  -  1  +  Aeb/a")) 

is  ^ 

-2_b/a.  Thus,  G(x|y)  <  e  /Y^-. 


NOTATIONS,  BILINEAR  FORMS,  SPACES,  AND  NORMS 

We  now  define  the  various  norms,  spaces,  and  bilinear  forms  used  throughout 

this  paper.  The  norms  introduced  here  are  analogous  to  those  defined  in  [2]. 
k 

The  space  H^(I),  k  =  0,1,...,  1  ^  p  <  00  is  the  usual  Sobolev  space  on  the 

interval  I  =  [0,1]  consisting  of  functions  with  k  derivatives  in  Lp(I).  On  this 
space  we  have  the  usual  norms  given  by 


luMHk(i) 

p 


r  z  f  iu<j)(,oi|’dx 

L  j=o  J  i 


Y.  ess.  sup.  |u 

J-0 


1  <  P  <  00 


°1  1 

The  space  H  (I)  denotes  the  subspace  of  H  (I)  of  functions  which  vanish  at  the 
P  P 

endpoints  of  I.  This  has  sense  because  H'*‘CC°(I).  Note  that  H°  =  L  . 

P  P  P 

Let  A  =  {0  =  xq  <  x^  <  ...  <  x^  =  1  },  where  N  =  N(A),  be  an  arbitrary  mesh 
on  the  interval  I  =  [0,1].  Let  h^  =  -  x^  1^  =  (x^  ^,x^.),  j  =  1,...,N, 

Pj  =  (hj  +  hj+1)/2,  j  =  1,  —  ,  N-l,  PN  =  hN,  and  h  =  max  Vu . 


We  seek  a  variational  setting  for  the  problem 


(2.7)  Lu  =  -eu"  +  a(x)u'  +  b(x)u  =  f  in  (0,1), 

u(0)  =  a, 

fu  =  gju'd)  +  82u(1)  =  8, 
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where  the  functions  a(x)  ,  b(x)  satisfv  assumption  Al,  and  flj ,  8-j  satisfy 
assumption  A2.  These  assumptions  Al  and  A2  will  be  assumed  to  hold  throughout 
the  entire  paper,  and  will  not  be  repeated  before  each  theorem. 


Let  L*  denote  the  formal  adjoint  operator  to  L,  i.e.. 


(2.8) 


2 

L*  =  -e-~  -  a(x)  4-  (b-a')(x). 
dx2  dx 


The  boundary  operator  adjoint  to  P  is  P*,  where  for  u  sufficiently  smooth 


(2.9) 


r.u-H  +  a<1,)U 


(1)  +  eu'(l),  if  Bx  +  0, 


u(l). 


if  B^  =  0. 


We  must  associate  a  bilinear  form  to  L,  and  describe  the  spaces  over  which  this 
form  is  defined. 

First,  we  define  the  space  H°  ,  1  _<  p  _<  °°,  to  be  the  completion  of 

p ,  A 


H,  =  (u€H  (I):  u(0)  =  0,  u(l)  =0  if  B,  =  0l, 

1  p  1 


with  respect  to  the  norm 


(2.10) 


IN!  „o 

hp,a 


1  Nl  1/p 

J*0  jupdx  +  I  |u(xj)  |P  ,  1  <_  p  <  °°, 

j  =  l 


II « I! 


where  N-^  = 


N  -  1,  if  Bi  =  0 


if  Bi  t  0. 


The  space  H°  .  can  be  easily  identified  with  L 

p ,  u  '  p 

u  =  (u,d^, . . .  ,djj^)  C  H°  A  =  Lp  ®  RN1,  and 


N, 

R  A,  that  is. 
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(2.1D  ||u||  o 

n,A 


[||uilL  (I)  +  Z  p  |d  |P]1/P  ,  1  <  p  <  co, 

P  3  =  1  .!  .1  - 


max  [  1 1  u 

j 


UU, 

°°  J 


In  consistencv  with  our  definition,  we  sav  u  C  H°  . ^  H^(I)  if  u  C  HP(I)  and 

p,A  p  p 


d.  =  u(x . )  for  i  =  1, . . . . ,Ni 


Note  that  the  norm  ||  •  |[^o  is  very  close  to  the  L  norm.  The  term 

P,A  P 


1 

l  P. |d. |P  is  the  trapezoid  quadrature  rule  for  the  function  |u(x)|p,  when  u(x) 
1  =  1  J  3 


is  continuous  and  d-j  =  u(x^).  Therefore,  for  any  continuous  function  u 

1  im 1 1  u  j |  o  =  21//p  | !  u  ]|  .  Because  of  the  boundary  layer  behavior  of  the 

Yr+o  p,A  p 

solutions,  an  Lp  (particularly  L^)  type  norm  is  appropriate  for  measuring  the 

errors.  The  quality  of  our  approximations  is  measured  in  the  H°  .  norm.  In 

p,A 

particular,  the  computational  results  and  adaptivity  presented  in  [35]  are  based 

on  the  H?  ,  norm. 

1»  A 

Let  us  also  define  H2  =  (v  C  H1 (I) :  v(0)  =  0,  v(l)  =0  if  8,  =  0,  and 

q,A  q  1 

v |  j  C  H^(I.) ,  j  =  1, . .  .  ,N },  for  1  £  q  £  00 .  We  will  equip  this  space  with  a 

j  9  3 

norm  to  be  defined  later. 

On  H  .  x  H  .,  where  —  +  —  =  1,  1  <  p  <  °°,  we  define  a  bilinear  form 

p,A  q  >  A  p  q  -  - 


(2.12)  B^Cu.v)  =  ..Ej  /  j  uL*vdx  -  ed.J(v'(x^))  +  dNr*(v), 

where  J(v'(xj))  =  v'(xj+0)  -  v' (x^-0)  for  1  <_  j  _<  N  -  1,  and  =  0  if  =  0. 
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The  limits  v'(x.+0)  are  well  defined  because  v|l.  £  H^(T.)  for  each  i. 

1-  J  q  .1 


Now  we 


will  furnish  the  space  H  .  with  the  norm 

q ,  A 


(2.13) 


sup 


j||  •  HI  ,  defined  by 
Ba(u,v) | 


uCHO  .  !l u II h° 


In  order  to  verify  that  l||  *  ||j  is  indeed  a  norm,  we  must  show  positive 

definiteness  -  linearity  and  the  triangle  inequalitv  are  evident.  To  prove 

positive  definiteness  we  will  use  the  identitv 

N  "  N-l 

(2.14)  v(y)  =  jJ1  G(x  |  y)  (L*v)  (x)dx  -  ,  J(v'  (x_1))G(xj  |y) 

+  G ( 1 | y)T*v. 

OO 

That  (2.14)  holds  for  v££  [0,1]  follows  from  the  properties  of  the  Green's 

2 

function.  By  a  density  argument,  and  using  the  fact  that  Hq  is  continuously 

embedded  in  and  hJ  together  with  the  fact  that  L*  is  a  continuous  mapping  from 

2  2 
Hq  into  Lq ,  it  follows  that  (2.14)  also  holds  for  each  v£  Again,  see  [34] 

for  a  more  detailed  proof. 

2 

For  v£  H  A  1  <  q  <  °°,  select  u  £  H°  .  ,  1/p  +  1/q  =  1,  where 
M  1  ^  *  —  -  U  p  j  A 

u  =  (u  ,  d. , . . . ,d„  )  and 
o  o’  1’  ’  nx 

u  !  =  sgn(L*v) | j  ,  for  j  =  1 . N, 

dj  =  -sgn( J(v' (x^ ) ) ,  for  j  =  1,...,  N-l, 

and  dN  =  sgn  T*(v)  if  Bj  4  0. 

Then 

N  N-l 

lBA(uo’v)l  =  j=]  ^ I j  |L*v|dx  +  .f1  E|J(v’(Xj))|  +  |T*(v) | . 
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If  B^(u0,v)  =  0,  then  L*vjjj  =  0  a . e . ,  J (v* (x^ ) )  =  0,  j  =  1,...,N-1  and 
T*(v)  *  0.  By  (2.14) ,  v  =  0.  Thus,  jjj  ’  |||  is  a  norm  on  ^  ,  for  1  _<  q  £ 

Let  us  conclude  this  section  by  introducing  another  norm.  For  any 
vCH^,  define  |J  •  ||H2  by 


N-l 


^  T  |L*vlqdx  +  ^  Eq|  J(v'  (xj))  |qP-~q 


(2.15) 


i  v !  L2 

q,e,A 


hN_qlr*(v) ^ 


1/q 


,  1  <  q  <  00  , 


maxf  max  ||  L*v|)  .  .  max  £ | J(v' (x.) ) | p, 1 , 

1<  j<N  l<j  <N-1  3  3 


|r*(v)h^],  q  = 


EQUALITY  OF  NORMS 

We  shall  now  prove  that  the  norms 


H 


and 


q  ,e,A 


are  equal. 


Lemma  2.6.  Let  vfH  .  then 

-  q,A 


(2.16) 


v  II  u2 

q,e,A 


1  <  q  < 


Proof:  That  ||  v||  „2  <  ll|v|||  follows  from  Holder’s  inequalitv.  For  the 

H  »  — 

q ,  E ,  A 

inequality  in  the  other  direction,  for  a  given  v,  u  =  uv  is  selected  such  that 


|vl||  =  sup„ 

uCH 


1  BACU, v) I  |Ba(uv,v) 


p,A 


I u I  I H° 

P,  A 


P,A 


M  L2 

q,E,A 
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If  1  <  q  <  0,1  this  is  done  bv  selecting  "  -  (u,  d.,....,d  ) such  that 

—  vi 

u  |  =  j  I.*v  |  q  ^sgn(L*v)lj  ,  1  <  j  N 

i  i 

d.  =  -  eq_1| J(v’ (x.)) !q"1Pjq"1sgn(J(v’ (x.))) ,  1  _<  j  <  N-l 

and 

dN  =  |r*(v)!q"1hj“qssn  (r*(v)>  if  ;-j  +  n. 


first  assume  |  jv|  |  2 


|L*v|  I  Let  r'  >  0  be  given  and  define 

Lco(Iy- 


En  =  (xO  7 :  |(L*v)(x)  !  >  !  i  L*v  | 


L  (I.) 
oo  J 


Then  m(Er))  >0,  where  m(A)  is  the  Lebesgue  measure  of  A.  Select  =  (u,0,...,0) 
such  that  u  =  (m(E  ))  ^sgn(L*v),  where  .<.  denotes  the  characteristic  function 


of  the  set  A.  Then  |  |  |v|  ]  |  _>  |  |v|  |„2 


-r:,  and  since  n  was  arbitrary, 


IIMII>IMIh2 

If  I  I  V I  I  2  =  e|j(v' (Xj)) Ipj1  or  ]  | v |  | H2  =  |r*(v)|hN1  the  selection  of  u^ 

°°,e,A 


is  obvious. 

0  2 

Note  that  Lemma  2.6  implies  that  is  1-1  proper  over  ^  X  ^  for 

1  _<  p  when  1/p  +  1/q  =  1. 


EMBEDDING  RESULTS 

The  following  lemma  which  is  a  slight  modification  of  a  result  from  [28]  js 
used  to  prove  the  embedding  result.  All  constants,  C,  ,  C2,...  appearing 
in  Lemma  2.7  and  Theorem  2.8  are  independent  of  p,  q,  v,  e  and  A. 

Lemma  2.7.  Let  G(x/v)  be  as  defined  in  (2.6).  Then 
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and 


9G  ( x/ v)  <  —  e  ^  +  C ,  f or  0  <  y <  x, 
3v  '  ~  ' 


3G  (x/v)  _<  —  e  ~(X  v)  /C  +  C9,  for  x  <  v  <  1. 

3  v  C 

Proof:  See  [28]  or  [34]. 

Theorem  2 . 8.  (Embedding  result).  If  v  £  fr  ^  ^  ,  then  v  £  Lm(I)n  H^(I)  with 


(2.17) 


|vllL  (I)  ^  cil  lvl  Ih2  i  ^IMIh2  .  •  1-q-"’ 

°°  l,e,A  q,c,A 


(2.18) 


l'"llL(i)  <V1/t,'1IWI„2  s 

q  q,f,A 


1  <  q  <  oo, 


where  C-^  and  C2  are  independent  of  v,  q,  £  and  A. 

Proof.  Inequality  (2.17)  follows  directly  from  (2.14),  (2.15)  and  Theorem  2.5. 
In  order  to  prove  (2.18),  let  G(x|v)  be  as  defined  in  (2.6).  Then,  for 


v  £  ,  we  have 

q,e,A 


v(v)  =  J'  G(x|v)  Kv(x)  dx  -  J~  G(x | v) c(x)v(x)dx 


-  £  c J(v' (x.))G(x. |v)  +  G(l,y)F*(v) , 

j=l  J  J 

where  ^(x)  C  Lq(I),  kJj  =  L*v|j.  ,  1  _<  j  _<  N,  and  c(x)  =  (b-a*)  (x) .  So, 

i  i 


(2.19) 


v’  (v) 


-f 


K  (x)dx 

dV  V 


/ 


7r~  (x j v) c (x) v(x)dx 

d  V 


-  £  e  J(v’(x.))  (x.|  y)  +  ~  (1  |y)r*(v) 
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=  w^y)  +  w0(v)  +  w_j(v)  +  w^(v) 


First  consider  w^  (v),  where 


:2-20'  =  /„ 


J  K  (x)dx  -  f  2%bi  K  (x)dx  +  /  atip 

n  9>’  v  Jo  3>’  V  J  y  ^ 


Kv(x)dx 


z1  (v)  +  (y) 


Using  Lemma  2.7,  we  have 


!  z 1 (y) 


V  l 

-'X 


ie-i<y-’0^|Kv(x)|dll  +  c2|  |kv|  !l_ 
1  . 

r  —  A  t— — 1  v  i'v'i  I  a  x-  -L.  r  II  v  II 


\  f  7  vC^r1)  iKv(x)  !  dx  +  C2  I  |Kv|  j^. 


where 


(x)  = 


e~ax,  if  x  <  0, 


0,  if  x  >  0. 


Extend  |K  |  by  0  to  all  of  R  .  Then  through  Young's  inequality,  we  have 


(2.2D  1 1  i:i  I  It.  (i,  -ci|Ii*(7)|Il  (r)!Ikv!!l  a,  +  c2N'<vllLin) 

q  1  q  1 


—  C  I  !  Kv  I  i  Lq  ( I ) ,  for  1  _<  q  _<  °°. 


Again,  using  Lemma  2.7,  we  have 


(2.22) 


'Z2^L  (I)  -  Cc  I  k  (I)’ 
q  i 

i 


x  g  q 

Also,  from  (2.19)  we  have  w  (v)  =  -  /  p—  (x I v) c (x) v(x)dx.  Using  a 

2  '  Jo  3v 

similar  agrument  to  the  one  used  to  bound  w^(v),  together  with  inequality  (2.17), 


we  have 


Young’s  inequality  states  that  for  1  <  q  _<  »,  if  sC  L  (R  ),  and  g  C  L^(R  ), 


then  h  *  s*g  exists  a.e.,  belongs  to  (R  ),  and 


L  ^  imil  m*ml  ■ 

q  q  i 


i 


(2.23) 


K'Il  (t)  a 

q  q  q  ,  C  ,  A 


Next , 


for  1  <  q  <"  ro. 


N-l  -  N-l 

,(>’)=-  £  eJ(v’(x.))  —■  (x.  jv)  =  £  z  (v), 

;  j=l  1  '  1  7=1  1 


f  |z  (v)|C,dv  eq  i  J(v' (x  ))  !q|  f  '  || 

J  0  1  '  (  J  0 


(xjy)l' 


y  i 


+J  !4~  (x.  |y)  |q|,  for  1  <  q  < 


By  Lemma  2.7, 


Therefore , 
(2.24) 


z.j  I  lL  (T)  _f  Cc1/q  1[c|  J(v’  (x.))  |  ]  ,  for  1  _<  q  _<  °°. 


!w3!Il  _<  V-  !  !z.f  (T  <  Cc1/o_1  £  cj  J(v'  (x.))  I  , 

q  t,  n  i-i 


for  1  <  q  <  °°. 


Finally , 


v4(v)  -  ^  ( i , v) r * ( v) . 


and  Lemma  2.7  Implies  that 
(2.25) 


„  I  -  C'1 /q  1  I r*\ 

1  w4  I.  (1)  -  “  '  •  ' 

q 


Expressions  (2.19),  (2 . 20)- (2 . 25) ,  and  the  triangle  inequality  implv  that 


|v  I  |'L  Cc1/M  ]  !  |  v  I  Ih2  ,  for  1  _<  q  _< 

q  q.f.A 
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exlstlnh:  of  sections 


Using  assumption  AT  and  Theorem  2.8,  it  follows  that  F  is  a  bounded  linear 

functional  on  .  for  1  q  -  .  Equation  2.  IT,  Lemma  2.6  and  a  simple 

verification  oj  (2.1)  implv  that  B.(u,v>  is  (1-1)  proper  on  H  x  11“ 

p,.:.  q. 

11  o  ’ 

for  1  <  p  <  » ,  where  —  +  —  =  1  .  For  1  <  p  ■  the  spaces  HL  ,  and  H“ 

-  P  q  P ,  2.  q, 

are  retlexivo.  In  this  case,  we  may  apply  Theorem  2.1  which  leads  to  the 
existence  of  a  unique  w  £  ,  solving 

(2.26)  F.(v,v)  =  F(v)  for  each  v  £  H“,, 1. 

'  q 

If  ;  *0  or  r^O  if  Bj-0  we  proceed  in  the  standard  way  writing  u=w+u_( 

where  u  is  a  piecewise  linear  function  on  ‘  which  satisfies  the  essential 
0 

boundarv  conditions,  and  w  G  HC  ,.  We  remark  that  since  u  is  bounded 

P»  • 

independently  of  •■,  w£  ,  and  !  j  w !  |  o  is  bounded  independently  of  r. 

'  9  -  ■  f *  , 

H  9  \ 

By  treating  the  boundary  conditions  in  this  way  we  can  restrict  the  theory 
to  the  case  of  homogeneous  essential  boundary  conditions  without  any  loss  of 
generality.  Therefore,  this  restriction  is  made  for  the  remainder  of  the 
theorems  in  this  paper.  These  restrictions  are  not  imposed  in  the  numerical 
examples  (see  [3T]). 

Another  assumption  which  will  be  used  later  is 
k+1 

A4  :  a ( x )  / 1  £  C  (T  .)  , 

b(x)/  C  Ck(I.), 

j  •' 

f  (x)/T  G  Ck(I .) ,  and 
°  'j  J 

a(x)  and  b(x)  are  independent  of  >  .  The  specific  value  of  k  will  depend  on 


the  finite  element  trial  space. 


CHAPTER  1 

L*-SPLINF  TEST  FUNCTION’S  AND  THE  INF- SUP  CONDITION 
In  order  to  obtain  a  finite  dimensional  approximation  to  the  solution  of 
(2.26),  the  finite  dimensional  spaces  S.  .  and  S0  .  must  he  specified.  To  obtain 
quasi-optimalitv,  these  spaces  must  have  the  propertv  that  the  inf-sup  constant, 

C 2,  of  Theorem  2.2  is  bounded  away  from  zero  independently  of  £  and  A.  UTien 
both  S^  ^  and  S0  ^  are  the  conventional  piecewise  polynomial  spaces,  this 

condition  is  violated. 

For  the  trial  space,  S  .,  we  take  the  space  S  =  (u  £  C°  D  H°  .  :u I T  is  a 

->- » A  r  p ,  A  1 1 . 

.1 

polvnomial  of  degree  £  r},  that  is.  the  usual  space  of  piecewise  polynomials  of 
degree  r.  For  the  test  space,  S?  first  consider  the  space  of  L*  splines: 

2 

(3.1)  S  =  (v  C  H  .  :  L,vv[  is  a  polynomial  of  degree  r-2  if  r  >  1, 

L*v| j  =0,  if  r  =  1 }  . 

i 

It  would  be  ideal  if  we  could  use  the  test  space  -  not  onlv  would  quasi¬ 
optimality  result  (Theorem  3.3),  but  also  the  nodal  errors  would  be  zero 
(Theorem  3.4).  Unfortunately,  since  the  functions  a(x)  and  b(x)  are  not  constant, 
it  is  in  general  impossible  to  determine  the  basis  functions  for  exactlv. 
However,  it  will  be  shown  (Theorem  3.3)  that  quasi-optimality  is  preserved  if 
basis  functions  can  be  found  which  are  sufficientlv  accurate  approximations  to 
the  "ideal"  basis  functions  of  ST  . 

(n) 

Suppose  we  can  find  basis  functions  {ijj.  .  }  -j  =  1,...,N;  (’  =  -l,...r  -  2, 

*  y  .1 


satisfying 
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-1 ,  i 


10. 


elsewhere. 


(n) 

;!  ,  .  (x. )  =  6  . 

-1,1  i  i ,  .1 


for  i , j  =  1 , . . . ,N. 


(3.2b) 


(n) 

L*'K  1  = 

,  i 


p^)‘ 


on  T  .  , 
.1 


el sewhere , 


,1  i 

_  max 

Let  ni  '  1-1, ...r-2 


v'f  <(x.)  =  0,  for  Z  =  0 . r-2 ,  i,i  = 


1,  .  (x)  !  |  ,  and  n  =  max(n.).  Denote  by  the  space 

*  ^  oo  i  -1  ^ 


spanned  bv  these  basis  functions.  Note  that  =  g 

(',) 


THE  INF-SUP  CONDITION  OVER  S  x  S, 

r  L 


We  will  now  show  that  the  inf-sup  constant  C?,  is  bounded  away  from  zero 
independently  of  e  and  A  when  S1  ^  =  S^,  and  S2  A  =  SL  ^  and  n  is  sufficientl> 


small  (independently  of  O* 

Before  proving  this  result,  we  need  to  define  some  additional  norms  over 
the  space  of  polynomials  of  degree  r.  These  norms,  as  well  as  the  basic  idea  of 
the  proof  of  the  inf-sup  condition,  are  taken  from  [2]. 


On 


the  space  of  polynomials  of  degree  r  over  the 
v  _  i 

b . (x-x)  ,  we  define  the  following  norms: 


interval 


[x , x+h] , 


f(x)  = 


DB 


m 


h 


p 


([ 


z  Kirh 

1=0 


p  .  lp+1 


1/p 


1  <  p  <  CO  , 


I !f ! !l  fi,s+h]  =l 


max 


ib.h1 i 


0  <  1  <  r 


r  \ffh  f(x)  (x-i)idxl  P 


1  1/p 


i=0 


p  i+p-  i 


.  1  <  p  <  °°, 


'L  [x,x+h]  = 
P 


max 
0<i  <r 


I fp^  (x-x)1dx[ 


i+1 


and 


/ 


(|f(X)  lp+|f(55+h)  |P)  h 


r-2  | p 
+  Z  jx 


+^f (x)  (x-x)  *dx  |  P 


i=0  ^pi+p-i 


1/p 


'l  [x.x+h] 
P 


1  <  p  <00, 


max  f I f (s) | , | f (x+h) | , 


/fhf(x)  (x-x)  1dx  | 


i+1 


,  i  =  0 , . . . r-2] , 

p  =  00  . 


Lemma  3.1.  | \ * | | 1  ,  [[■[!"  ,  and  | | * | | 

a  constant  C  =  C(r),  independent  of  h,  f,  and  x,  such  that 


are  norms  over  S^,  and  there  exists 


(3-3a^  C  f x , x+’n  1  -  [x.X+h]  [x,x+h]  , 


(3.3b) 


P 


<  cl  I f M 

L  [x,x+h]  —  1  1  [x,x+h] 

P  P 


and 
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(3.3c) 


[x,x+h]  1  | | f | | L  [x,x+h]  X  Ci ] f j  |  L  [x ,x+h] , 
p  P  p 


The  proof  of  this  result  can  be  found  in  [2], 

We  will  now  prove  the  main  result  of  this  section.  All  constants  C,  C ^ ,  C9 , 
Dq,  Dj  appearing  in  the  theorem  or  the  proof  are  independent  of  c ,  .1,  p,  v,  and  n. 

Theorem  3.2.  Let  n  be  sufficiently  small  and  —  +  —  =  1. 

-  p  q 


inf^  .  sup  | B .  (u,v)  |>  Do(l-D1n) 

v£S,  uC  S 

L  r 


i v 1 1 h2  = ;  iHih°  = 1 

P  t  ‘  Pn- 


for  1  <  p  <  “. 


Proof .  Let  v£S^'  be  given.  In  terms  of  the  basis  functions 


N  N  r-2 

vCx)=L  v(xj)v-i,j(x)  +  £  £  b£,ihi’^,i(x)- 

i=i  j=i  1= o 


Then,  by  (3.2a, 

b) 

r-2  _  r-2 

(3.4) 

L*v|x 

.1 

=  2-  b .  .(x-x.  .  )*  +  b. 

i-0  3-1  i=0  1 

+  (v(xi-l)n-l,.i-l  +  V(Xj)n-l, 

Hence , 

(3.5) 

1  '  L*v  | 

1 

'L  (I.) 
q  i 

r-2 

r-2 

- 

|:L 

i=n 

bi,i(x-xj-i)1l!i.  (I.)  +  !!.^n  bi 
q  J  1=0 

+  :|v(x1_1)n_1>1_1  +  v(x.)n_Ml|^(I_v  for  1  :  q  : 
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Since , 


I  I  r-T 

lit  _ 


Z  h •  •  (x_x .  i ) 1 1 ' 

|i-0  1  lkq(I,) 


£  |b4  ,|Vi+1 

i=0  ''  1 

(r-l)nax|bi  .  |h*, 


,  1  <  q  <  oo, 


Z  b,  .h!n.  . 

,_n  i-J  .1  1.1 


L  (I.) 

q  i 


•i.  L  |b(  . t^h?1+ 
3  [i-o  1 

i1 .  (r-l)max  lb .  .  |h^  , 

1  ,  i,J  2 


,  1  <  q  <  oo. 


|v(xi-i)r,-i,j-i  +  v(x3),1-i,kli.q(i1)-  2r,j*  >v'  lLoi;(iJ)hi/q’ 

1  <  q  <  oo, 


it  follows  from  (3.5)  that 


(3.6a) 


Zll^vllJ  <  2q(l+n)q  E  £  |b  |V1+1 

i-1  q'  j1  j-1  i-0  ,;l  1 


+  «V||v||J  , 


1  <  q  <  oo  , 


(3.6b) 


max  L*v 


L  <  (l+n)  (r-l)  max  max  |b  .  |h! 

oo^V  l<j<N  0<i<  r-2  3 


+  2n 1 1' 


From  (3.6a,b)  and  the  embedding  result  (2.17),  it  follows  that  if  n  is 
sufficiently  small 


2  7 


0.7) 


E  E  Ik  JVi+1  +  E  ct’|j<v'(»,»qn,-<l 

i=i  i=n  1-1  5  i=i  ’  1 


^  ,  A 


+  !  r*(v)  iqhjJ.  q 


l/q 


,  l  <  q 


C  max[|hL<1|hIuc]J(v’(xi.))|p'1,  |r*<v)‘  l^,1]  , 


V 


where  L  and  J  are  the  indices  such  that 


(3.8) 


bT  , I h*2  =  max  max  |h,  .  |h., 

L'J  J  l<j<s  0  <r-2  *•’  ■’ 


and  K  is  the  index  such  that 


(3.9) 


£ | J(v’ (x  ))  j p  1  =  max  e| J(v’ (x.)) Id.1 

K  K  l<i<N-l  1  -1 


.  (n)  . 


With  v  £  Sj^  given,  select  uv  =  ^  +  <f7  in  the  following  wav.  Select 


*1  Sr  such  that 


(^(x^)  =0,  .1  =  0, . N,  and  if  r  >  1, 


(3.10) 


f  (x-x  1)1(t>1 
I .  J 


dx  =  ih,  .  |q  1hql+1  sgn  h  .  , 


1 ,1  i 


i  <  i 


for  i  =  0,....,r-2,  j  =  1,...,N,  1  ^  q 


Select  d>„  £  S  such  that 
2  r 


(3.11) 


i7(x.)  * 


q-i 


y  e  ;J(v’(x.))|q'1p!'q  sgn(J(v'(x.))), 


1  <  j  <  N, 
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and  if  r  ^  1 

f 

That  and 
Recall, 

(3.12) 

Bv  (3.4),  (3 


Therefore , 
(3.13) 


^(x^)  =  |r*(v)  Iq  sqn (T*v) 


1  <  q  <  »■■, 


(x-x.  ^)  $  f  =  0,  for  i  =  n,....,r-2  and  j  =  1,....,N. 


<J> 0  are  uniquely  determined  follows  from  Lemma  3.1. 


N-l 


Ba(u,v) 


fj  u(L*v)dx  -  ej(v'  (x.))d. 


+  d  r*(v). 
N 


10)  and  (3.11),  we  have 


BA(fj  +  <j!0,v) 


N  r-2 

-  z  z  i\  ^ 


q  ,  qi+1 


N-l 


;))|q 


£  eq  |  J(v '  (x .  ’ 

J-1  i=0  *’■'  •’  1  =  1  J 

N  /r-2  \ 

+  |r*<v)  |"  +  £  /j.ti  h./fij  <*i  + 

N 

+  Z  f  -i  +  v(xj)0  1  +  1  _5  q 


i>2)dx 


Ba(*i  +  4>2»v) 


N  r-2 


N-l 


,  E  E'  h.  < !q  t,?1+1  +  E  c',U<v(x.))|<’ 

-  1-1  i =0  1,1  1  J-1  1  ’ 

+  ir*(v)  iV~q 


“  n 


N  r-2 


E  E  hj  Jq  i>11+1 

j-1  i-0  1,1  1 


1/q 


+  2 


IVHl„|  I  1*1  +  ^2>!  o 


1  <  q  < 


co 


From  the  definitions  of  £  and  ((3.10),  (3.11)),  (3.3c),  and  (2.11) 


we  have 


(.3.14) 


N  r-2 


1 ! :  +  <  cp  £  £  !ii.  .I"  i,?,+I 

•  1  2'  o  -  3  =  1  1=0  i,.l  1 


P,A 


+  Y,  cq I  J(v’ (x.)) |qol'q  +  !r*(v)!q  hA  q  , 

1=1  1  1 


1  <  p  <  oo 


Inequalities  (3.13),  (3.14)  and  (3.71  vield 


B.o,  +  -:2.v) 

Tit i  +  «,H  „ 

-  H 

r*A 


c,HvP  2  -  c  nllvll 

p  ,e ,  A 


>  Dj  (l-D.7Ti)  !  [  V I  I  2  .  1  <  q  <  00 , 

q,e,A 

with  the  last  inequality  following  from  the  embedding  result  (2.17). 


Next,  consider  the  case  when  p  =  1  and  q  =  °°.  In  this  case,  we  modify  the 

if  i  4  J, 


definition  of  <i^  in  (3.10)  such  that 


f  (X“X.  ,  3 14>1dx  = 

T  1  1 


I  . 

.1 


o. 


5i,IhJ  SSn(bL,J)’  lf  j  =  J’ 


where  L  and  J  are  defined  in  (3.8)  and  6  is  the  Kronocker  delta.  The 

i*J 

definition  of  cp  in  (3.11)  is  also  modified  such  that  <f>„(x.)  =  6.  0  "*sgn  J(v'(xK)), 

l  L  J  .1  »K.  K  N 

for  1  _<  i  N-l,  and  (^(x^)  =  sgn  T*(v),  where  K  is  defined  in  (3.9) 


Bv 


(3.3c)  I  !£>.,+<£  |  I  <  c.  From  the  modified  definitions  of  and  4>-i»  (3.12) 
i  i  j|0  —  1  l 

1 ,  A 

and  (3.7)  it  follows  that 

B  A(<P,  +  <f>9,v) 

- - - -  <  D1(l-D2n)  IM  |  2 

I  1*1  +  *2II  0  “ 

hi,a 


as  desired. 
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When  v  =  «*  and  q  =  1  the  formulas  for  tj  and  given  in  (3.10)  and  (3.11)  are 


still  valid,  and  from  (3.3c)  we  have  jj*  4*  ]  j  q 

1  —  H 


K'l  +  «2I  lL  i  C.  The 


desired  result  then  follows  from  this  fact,  and  inequalities  (3.13),  (3.7)  and 
(2.17) 

Theorems  3.2  and  2.2  yield  the  following  result. 

Theorem  3.3.  Let  u  C  H°  ,  he  the  solution  to  B.(u.v)  =  F(v),Y  v  f  , 

— - - —  A  q  r  ,A 

1  <  p  <  °°,  —  +  —  ~  1.  Furthermore,  assume  u  C  Then,  for  r|  sufficiently 

p  q  °°,A 

small,  there  exists  a  unique  solution  u^  £  to  the  problem  B^(u^,v)  = 

F(v),V  v  C  and  |  | u  —  il  ]  |„0  <  C  inf  |  |u  -  w|  |  0  ,  1  <  p  <  “=,  with  C 

t.  L  H  o  H  . 


P ,  A  vi 


independent  of  u,  p,  e,  and  A. 


The  next  theorem  shows  that  if  q  =  0,  then  (u  -  u  )(x.)  =  0,  for  i  =  1,...,N. 

L  J 

Theorem  3.4.  Let  u  £  H®  .  ,  1  <  p  <  »,  and  u„  £  S,.  C  H®  ,  ,  with 
-  p,A  —  —  L  1  p ,  A 


(3.15) 


BA(Vv)  =  BA(u,v),  V 


Then  u  (x. )  -  d  =  0,  i  =  1, . . . ,N  ,  where  u  =  (u,d. . . dN  ) 


(3.16) 


O 

Let  1  <  q  <  00 .  Bv  Theorem  2.1  there  exists  v.  (H  .  such  that 

1  q,e,A 

B,(u,v.)  =  d.  V  u  £  H°  .  ,  -  +  -  =  1. 

A  i  i  p,A  p  q 


From  (3.16)  and  the  definition  of  B^  (3.12),  it  follows  that  L*(v)  =  0  on  every  £, 

J(v'(x  .))  5  .  for  i ,  j  =  1, - -  N-l , 

J  J  ^  1  5.1 


r«(v  )  =  a 

v  V  i,N 


i  =  1 , .  . . ,N. 


Because  v^(x)  is  continuous,  and  L*(v)  =  0,  on  (0,x^)  and  (x^,l),  is  the 

Green's  function  at  x  =  x, .  Thus,  v  £  ,  and  also  v.  £  This  implies 

that  (3.16)  holds  for  1  <  q  <  °°,  and 
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4 


Vv^V  =  n  =  ul(xi)  -  di’ 

which  finishes  the  proof. 

Theorem  3.4  is  a  restatement  of  the  well  known  fact  that  when  the  Green's 
function  at  x  =  x^,  x^  C  A,  belong  to  the  test  space,  then  the  error  at  the 
nodal  points  is  zero. 

As  pointed  out  earlier,  the  space  is  in  general  unobtainable.  However, 

in  light  of  Theorem  3.3  it  will  be  satisfactory  if  we  can  generate  the  basis 
functions  for  the  space  \  (3.2a,b),  provided  p  is  sufficiently  small.  This 

will  be  done  in  the  following  chapter. 
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UNCLASSIFIED 


CHAPTER  4 


EXPLICIT  REPRESENTATION  OF  THE  TEST  SPACE-S 


(n) 


In  this  chapter,  we  will  determine  the  basis  functions  (3.2a,b)  for  S 


explicitly,  with  the  condition  that  n . 


I  In 


(n) 

L 

k 


max  |  •  I  lL  (I  )  1  Ch  The 

i=-l . r-2  1 


constant  C  will  depend  on  l'min  a ( x )  ,  and  the  local  smoothness  of  a(x)  and  b(x) 

1  . 

! 

on  Ij,  but  will  be  independent  of  h  |  and  e. 

In  order  to  determine  the  basis  functions  for  in  Ij  ,  we  first  rescale  the 

interval  Ij  to  I  =  [0,1],  and  then  drop  the  index  j  for  simplicity.  A  "tilda" 
will  be  used  to  denote  this  rescaling,  for  example,  if  g(x)  is  defined  for 
x  £  Ij,  5 (v )  is  the  function  defined  for  v  £  I  such  that  g(v)  =  g(yhj  +  x.^). 
Recall  that  c(x)  =  b(x)  -  a'(x). 

After  rescaling,  our  goal  is  to  seek  approximations  to  the  solutions  of 


(4.1a) 


L*v 


~  v'  +  c/(y)v  =  in  I, 


v(0)  =  1,  v(l)  =  0, 


(4.1b) 


L*v  =0,  in  I , 


and 

(4.1c) 


v(0)  =  0,  v(l)  =  1, 


L*v  -  y  .  in  I,  Z  =  0,...,r-2 


v(0)  =  v ( 1 )  =  0. 
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UNCLASSIFIED 


Because  of  our  requirement  on  the  space  in  (3.2a,b),  we  must  estimate 

ML  i,.i 


Notice  that  the  rescaled  operator  ha?  the  coefficient 


in  front  of  the 


highest  order  derivative.  If  a  standard  asvmntotic  expansion  in  powders  of  €  is 

used  (see  e.g.  [10],  [15],  [36]),  then  the  errors  ^an  be  shown  to  be  no  better 
£  k 

than  0((v)  ).  This  is  undesirable  because  one  of  our  goals  is  to  use1  adaptive 
mesh  refinements  which  will  quickly  lead  to  intervals  having  size  h,  with  h  <  fc  . 

Although  the  asymptotic  expansion  in  c  will  fail  as  h  becomes  small,  it  will 
be  shown  that  an  asymptotic  expansion  in  h  will  produce  errors  of  size  O(h^) 
independently  of  e.  In  order  to  validate  our  asvmntotic  expansion  the  following 
lemmas  are  needed. 

Lemma  4.1.  Let  v  solve 


r  v"  +  iLiili  v>  - 


g  in  1 , 


v(0)  =  v ( 1 )  =  n, 


with  a(v)  >  a  >  0  and  g  bounded. 


!  | v 1  I  <  min  (h/a,  2h  /e)|jg|> 


Proof.  From  the  maximum  principle  ([32]  Chapter,  1,  Theorem  11),  it  follows  that 
! v ( v )  <  | w i ( v ) |  where  w,(y)  =  |!g!j  h(l-v)/a,  and  |v(v)|  <  |w?(y)|  where 

L'cr  i- 

Wo(v)  =  h2(e-ev)  |  |g|  I,  /e. 


Lemma  4.2.  Let  w  solve 


w  =  g 


w(0)  =  0, 
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where  aG  '■  0  is  constant  and  g  is  bounded.  Then  |  |wj  j  |  <  h(l-e  a°',C)|  |g!  |T  / 


Proof .  This  result  follows  immediately  from  the  identity 

•>  -v 


w(v)  = 


-  -  T  f  • 


-u,,h(s-v)/c 


o 


P(s)ds* 


0 


(n) 


We  seek  an  approximation  v  to  the  function  v,  where  v  is  the  solution  to 
(,1.1a),  (4.1b)  or  (4.1c).  Assuming  sufficient  smoothness  of  the  coefficients 
a(v)  and  c(v),  we  may  expand  bv  the  Tavlor  series  around  v  =  0  to  obtain 


(4.2) 

where 

(4.3a) 


(4. 3b) 


L*v  =  -  -~r  v"  -  — ^ -  v’  +  c(v)v 
h2  h 


-  V"  -  A(2>  v- 
2  h 


+  h[-  v’  +i^-  v] 

n  n 


v. +vcy<0)v, 


k  a<k>«  > 


,  ,  kr  Z_ _ a_  ,  y 

h  k!  h  (k-1) ! 


k-1  c(k_1)(r  ) 


—  v]  , 


where  f.  ,  C  C  (0.1) 

3  C 


Therefore,  we  can  write 


L*  = 


h,  +  hLi  +  ' "  +  •’’"Vi  +  hV 


c  d-  1(0)  d  , 


J°  h2  dv2  h  d>7 


L  =  i  r  _  £  a(-1)(0)  A  .  .^:L  c(i-1)(0)l 

1  h  1  j!  ' ’  dv  (i-1) !  (U;|  ’ 


i  =  l,...,  l-i  , 
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( A . 3c) 


£  t  - 


f 

y 

Z\ 


( o 


We  seek  in  the  form 

..(n) 


=  v  +  liv,  +  h“v„  + 
o  1  2 


(4. A) 

From  (4.2)  and  (4.4)  it  follows  that 


,  .  k-1 
+  h  v 


k-1 


?.  =  1 


>•5) 


l*v(n)  -  L  v 
o  o 


+  hi  L0vx  +  1.^1 


+  h  fl.  v  +  L. v.  +  L,v  l 
o  2  11  2  o 


k-1 


+  h  f  I,  v  +  L.  v  0  +  •  ■  •  1.  v  ] 
o  k-1  1  k-2  k-1  o' 


+  h  [ K1vk_I  +  Vo:  +  +  VJ 


Based  on  this  formula,  the  functions  v  ,  v  ,  ...  should  be  defined  recursively 


°  1 

follows : 

(4.6a) 

L  v  =0, 

i f  a'  ^roximat ine 

o  o 

the  solution  to 

c 

II 

V — ' 

G 

> 

II 

C 

c 

> 

(4.1a)  , 

Ah) 

C 

II 

0 

> 

c 

if  approximating 

the  solution  to 

VQ(0)  =  0,  vp(l)  =  1 , 

(4.1b), 

as 
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(4.6c)  L  v  =  v\  l  =  0 . r-2, 

OP 

v  (0)  =  v  (1)  =  0,  if  approximating 

o  o 

the  solution  to 
(4.1c), 


i 

I 


L  v .  =  -  /,  L  .  v .  . , 

°  1  ft  J  ^ 


(4.7) 


Vj (0)  =  v.. (1)  =  0, 


i  =  1, . . . ,k-l . 


From  (4.7)  and  (4.5)  it  follows  that 


(4.8) 


L*v(fl)  =  L  v  +  hk  [  £  R  v  .] 
o  o  i  k-iJ 

i=l 


This  leads  us  to  the  following  theorem  which  is  the  main  result  in  this  chapte 

.  ~  k-1 

Theorem  4.3.  Assume  that  a£  CK(I),  and  b  C  C  (I).  Let  be  defined  by 

either  (4.6a),  (4.6b)  or  (4.6c),  and  v . ,  i  =  l,...k-l,  defined  recursively  by 

(4.7).  Let  v(f!  -  v  +  hv.  +  .  .  .  +  ^v.  .  .  Then 

o  i  k-1 

||L*v<rl,||L„<  eh1'’1 

if  v  is  defined  bv  (4.6a)  or  (4.6b),  and 
o 

! |L*v(n)  -  L  v  |  L  <  Chk 
■  1  o  o  1  1 L  — 

CO 

if  v  is  defined  bv  (4.6c),  and  in  each  case  C  is  independent  of  e  and  h. 
o 

Proof.  First  assume  v  is  defined  bv  (4.6a).  Because  of  equation  (4.8)  it  must 
k  ° 

be  shown  that  h| j  ]T)  R.v  . ]|„  <  C.  We  prove  this  by  induction  on  k. 

•  _ ,  IK-1  L 
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First,  take  k=l.  From  (^.3c)  (with  £.=  1 ) ,  we  have  hR.v  =  -va'(£  )v*  +  c(£  )v  . 

1  o'  a  o  c  o 

From  (4.6a)  and  the  maximum  principle,  it  follows  that 


(4.0) 


0  <  v  (v)  <  1 
—  o  — 


Let  w(v)  =  vv'  -  (v  -  1).  Then 
•  o  o 


_  w-  .  w  =  ±(0)  (  _  1} 

.2  h  h  ',vo  1 '  ’ 

h 


and  w(0)  =  0.  By  Lemma  4.2  and  0.9)  we  obtain  |  |w|  <  1.  Hence, 


vv 


o  1 L 


[w  +  (vQ-l) | I  <  2,  which  implies  that 


! hRivo 1 1  1  +  ! Ic| lL  i  C 


with  C  independent  of  e  and  h.  Thus,  our  assertion  is  true  for  k  =  1. 
Next,  assume  that 


I v± |  | L  £  C,  i  =  0,1, . . . ,k-l  and 


(4.10) 


I  yv! | ! 


i  ML 


<  C,  i  =  0,1 . k-1. 


This  is  actually  our  induction  assumption  because  from  (4.10)  it  follows  that 
|hR^_^v^||  _<  C,  for  i  =  l,...,k-l.  We  must  show 


and 


!vkML 


|yv’llL  <  c- 


From  the  induction  assumption,  (4.10),  the  definition  of  v^»  (4.7)  and  (4.3b) 


we  have 
(4.11) 


lLo\HL 


k 

£  LiVi 

'i=l 


i  £  * 


and  so  bv  Lemma  4 . 1 
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(A. 12) 


ic- 


As  before,  let  w  =  yv^ 


v,  .  Then 
k 

£  ,  3(0) 

-  -  w  —  — - - —  w 

,2  h 


t  .  a(0) 
vL  v,  +  — - —  v,  , 
ok  h  k 


and 


w(0)  =  0. 


Inequalities  (A. 11),  (A. 12)  and  Lemma  A. 2  imply  that  | |w| £  C.  Therefore, 

00 

| |yv^| £  C  which,  with  (A. 12)  implies  that  | |hR^v^| £  C.  This,  together 

OO  00 

with  the  induction  assumption,  proves  the  result  when  v^  is  defined  by  (A. 6a). 
When  vq  is  defined  by  (A. 6b)  the  proof  is  almost  identical. 

When  v  is  defined  bv  (A. 6c)  we  use  Lemma  A.l  and  deduce  that  I  |v  II  <  Ch. 
o  -  1  1  o 1  1 L  — 

00 

If  we  set  w  =  yv'  -  v  and  use  Lemma  A.  2  it  will  follow  that  |  |w'  ]  |,  <  Ch. 

o  o  1  1  •  o 1  1 L_  — 

00 

Induction  on  the  assertions 

I | vi I lL  <  Ch, 

00 

and 


i  L  — 

OO 


yields  the  desired  result. 


This  theorem  proves  that  it  is  possible  to  choose  the  basis  function  f in 

i.J 

such  a  wav  that 


|l<j  .  .  Iln^jll^  ,  iohj'1 


where  C  =  C(a|^  ,  b | ^  ,  k)  is  independent  of  £,  and  h..  Recall  that  by  Theorem 

\i  i  3 

3.3  a  small  value  of  n  will  guarantee  a  quasi-optimal  finite  element  solution. 
Since  we  are  approximating  the  exact  solution  with  a  piecewise  polynomial  of 
degree  r,  the  accuracy  of  the  method  will  not  increase  in  order  as  we  increase 
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Therefore,  it  is  sufficient  to  take  k  =  2,  in  which  case 


the  order  of  n. 

n  =  0(h). 

/  _  \ 

These  basis  functions  T.1.  are  easily  derivable  and  their  explicit  formulas 

i,l 

—  V 

are  given  in  [34].  In  general,  they  have  the  form  P^(y)  +  P?(v)e  '  ,  where  P^ 
and  P2  are  polynomials,  and  X  =  aji/e  is  °ften  referred  to  as  the  local  or  cell 
Peclet  number.  For  large  values  of  X  these  basis  functions  themselves  exhibit 
boundary  laver  behavior.  When  the  value  of  X  is  small  these  functions  are  close 
to  polynomials. 


However,  the  nodal  errors  will  decrease  as  n  decreases. 
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CHAPTER  5 

PROJECTION  ONTO  UPWINDED  POLYNOMIALS 

When  the  space  ^is  used  as  the  test  space  for  the  finite  element  method 

Theorems  3.3  and  4.3  yield  quasi-ontimality  for  the  approximate  solution  in  Sr> 

However,  since  the  basis  functions  of  S^^are  of  the  form  P^(v)  +  P2(v)  e 
aGh 

with  X  =  — —  ,  the  bilinear  form  requires  the  integration  of  functions  with 
boundary  layers  and  smooth  functions  as  well.  Unless  a  special  quadrature  rule  is 
used,  which  integrate  terms  of  the  form  vm  e-^'’  exactly,  large  quadrature  errors 
will  result  whenever  X  is  large.  Standard  quadrature  is  also  needed  for  the  smooth 
terms. 

In  order  to  avoid  this  inconvenience  we  propose  to  project  these  "exponentially 

upwinded"  basis  functions  of  Sj/^onto  a  space  of  polvnomially  upwinded  functions. 

These  projected  basis  functions  will  have  the  form 

M 

(5.1)  iMx)  =  $  (x)  +  ot  g  (x)  for  xCl. 

1  J  i=1  11  J 

where  <f.(x)  is  the  standard  piecewise  linear  "hat"  function  and  =  a^(e,l  '  are 
the  upwind  parameters.  This  is  a  direct  generalization  of  the  commonly  used 
a-quadratic  upwinding  in  which  case  M=1  and  g^(x)  is  quadratic  on  I  j  .  In 
particular,  it  was  shown  in  [  3  ]  that  when  b(x)  H  0,  and  a(x)  is  constar  the 


value  a,  computed  from  projecting  the  space  is  identical  to  the  so  called 
"optimal"  value  of  a  presented  in  ([8],  [9],  [17]- [21],  [30],  and  [37]),  which 


yields  the  exact  nodal  solutions  when  f (x)  is  constant. 


Let  .  be  as  defined  in  (3.2a,b)  and  as  constructed  in  (4.6a,b,c). 
J 

Define  s(^)  by 
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(5.2a) 


.00 


S^'7  =  Span^X^  .J.  j  =  l,  ....  N;  *— 1,  r-2. 


where 

(5.2b) 


(n)  k  (M) 


with  denoting  the  projection  operator  onto  the  first  k  Legendre  polynomials 


$Q,....,<|>b  on  each  interval  I...  Also,  we  take  the  convention 


1 


(5.2c) 


.(h) 


Xj(x)  f  X£,i(t)  +  ^,j(xj-l)’ 


Xj-1 


and  hence 


xM(xj-i}  ‘ 


and 


*M(V  ’  *M(V' 


Prom  the  results  of  Chapter  4  these  upwind  coefficients  t.  =  (^’ ,  4>.)j  can  be 

i  1  ,  i 

computed  exactly,  since  onlv  integrations  of  the  form  f  ymeAXdv,  or  integrations 

J0 


with  polynomial  integrands  are  required.  Once  these  coefficients  ,  i  =  l,...,k 


are  computed  on  each  interval  I.,  all  integrations  remaining  will  be  of  the  form 


J0  P(y)g(v)dy,  where  g(v)  is  smooth  independently  of  e. 


Let  vT  C  S  n  .  Then  P,  (v')  =  v  is  a  piecewise  polynomial  of  degree  k  with 

L  L  K  L  cl 

the  property  that 


(5.3) 


1  i 

(g,v--v  )  =  0,  for  each  g(x)  C  Sk, 


a  piecewise  polynomial  of  degree  <_  k.  Note,  that  since  va(0)  =  v^(0)  =  0,  it 
follows  from  (5.3)  that 


(5.4) 


vL(Xj)  =  va(Xj),  j  =  0, - ,N. 
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Before  proving  the  main  result  of  this  section,  we  first  prove  two  lemmas. 

.  (n) 


Lemma  5.1.  Let  vT  £  S.  .  Let  P,  (v* )  =  v'  where  P,  is  the  local  L„  proieetion 

L  L  K  L  t  K  Z 

operator  onto  piecewise  polynomials  of  degree  k.  Then 

l|vL  *  \I|L1(I>  i  |  !v; 

k 


V 1  1 V  lL1(I)  ’ 


where  C.  =  I  /2i+l . 
k  i=0 


Proof.  Let  4>^,  i  =  0,...,k  denote  the  Legendre  polynomials  of  degree  i  on  the 


interval  I..  Then 

j 


k  _ 1 


a  1 .  .  n 

.1  i=0 


=  v 


<t>.  , 


(*.,  40 


l  i 


where  denotes  the  L9  inner  product  on  I. 

J  -* 

Therefore , 


,  k  I  k.  I  L  fT  \  I  I  It  (1  \ 

7al  lLl(I  )  1  NVlI  lLl(I.)  Yj  - — 5 - 

J  _  n 


i=0 


If  we  normalize  the  Legendre  polynomials  in  such  a  way  that  =  1,  then 


H^Hl  a.)  =  and  (ij  =  2i+r 

OO  J  7  2  1 


and 


hence,  I  l*tl  1,.^)  £  h*/2|  l^l  !Lj(Ij)  .  Thus, 


k 

£ 

i=0 


I  ^i'  k  (I.)  I  ^i  I  k.  (I.) 


!  16. 


I  1'i'il  l^Clj) 


<  /2i+r 

i=0 


=  C, 


So 


r'  ! 


'a1  lL1(Ii)  -  "k1  |VL>  'Ljd  ) 


4  3 

UNCLASSIFIED 


1 


K. 


M 


It 


UNCLASSIFIED 


and  consequently. 


'V^d)  — ■  Ck  '  VI.  ' L^  (I) 


The  lemma  now  follows  from  the  triangle  inequality. 

Lemma  5.2.  Let  w  f  S„.  For  each  v,  £  ,  let  v  £  be  such  that 

-  r  L  L  a  n 

v'  =  P.(v'),  with  k  _>  r.  Assume  aj  £  C^+^(I.),  and  b|,  £  C^d*)  for  j 

Oi  K  L  X,  J  1  i  J 

j  J 


=  1 , .  . .  ,  N . 


Then 


sup  i  C  hri_r(  1  la<k+1'r)  I  l„r„  , 

vLcsL<">  i  J 


+  I  lb 


(k-r) 


for  1  <  q  < 


Proof .  Since  v^(x^)  =  v^(Xj)  for  j  =  0,...,N,  we  can  rewrite  the  bilinear  form  as 


BA(w-Vva}  = 


N 

£  / 

i=i  i. 


£w' -aw- 


x 

f 


(b-a' )wdt 


J-l 


(v'-v')dx 
L  u 


Because  of  (5.3)  and  the  fact  that  r  <  k,  we  have 


B*(w,v  -v  )  = 
A  L  a' 


N  r 

Z  / 

.1=1  I. 


-aw  - 


X 

/ 


(b-a')wdt  -  p(x) 


>1 


(v’-v')dx, 
L  a 


where  p(x)  is  any  piecewise  polynomial  of  degree  <k.  Consequently, 


(5.5) 


BA^W,VL~Vot) 


<  max[| |aw-p1(x) | |  ,  ,  +  /  (b-a ' ) wdt-p, (x)  (  J 

J  00  .i'  Jx.  ,  "  00  .V 


'3-1 
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Let  p^(x)  be  the  k  interpolant  of  aw  on  each  1^. 


That  is. 


Pl(x.  j+ihj /k)  =  (aw) (x  j+ih.  /k) 


for  i=0,...,k,  and  j=l,...,N. 


Then, 


(5.6) 


|av  -  Pl  (*)||La>(I  )  1  T(kTlTT 


r  ||(aw)(k+1)||Loo(I.) 


From  Leibnitz*  rule,  we  have 


(k+1) 


E  <  T> 

i=0 


k+1.  (k+l-i)  (i) 

,•  )  a  w 


V*  ,k+l.  (k+l-i)  (i),  since  k  >  r 

/  >  l  i  )  a  w  — 

i=0 


Hence, 


I  |  (aw) 


^  Il^Ij)  -  (r+1)(k+1)!!  Ia'  )  I  lH^(ij)  l^ll^dj), 


<  C(k+l)!||a(k+1”r)||  r(  sh  r| 

rr.  '  J 


'L  (I.), 

°°  J 


the  last  inequality  following  from  the  inverse  theorem.  From  this  and  (5.6)  it 
follows  that 


(5.7) 


law  -  Pl(x)| |L  (  }  <  rhk+l-r||a(k+1-r)||Hr (]>)l|wl|1  (T 

no  j  co  l  cr  1 


th  /*A 

Similarlv,  if  p»(x)  is  the  k  interpolant  of  J  (b-a’)w  on  I,, 

X 1-1  1 


(5.8) 


/  (b-a')wdt  -  P2(x)||  ( T  ( 
Xi-1  J 
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a 


h 


h 


N 


;  Chf1-r||(b-a')<k-r)|lHr(I  , 


‘L  (I.)’ 

00  y 


From  (5.5),  (5.7),  (5.8),  Lemma  5.1,  and  the  embedding  result  (2.18),  it 
tollows  that 


B,(w,v  -v '  1  „  uTi-r ,  |  | 

A  L  a  <  C  max  h.  { |  I  a  j  |  r 


sup 

V,  £  s,(n)  1  |vL  I  lH2 

L  L  q,  e,  A 


k+l-r , | |  k+l-r) i 
j  (||a  1  1 H 1  ( I j  ) 


+  ||b(k'r) 

We  are  now  ready  to  prove  the  main  result  of  this  chapter 


1 H  (I.) 
”  2 


MV' 


Theorem  5.3.  Suppose  that  assumptions  A1  -  A4  hold,  and  that  S  is  defined 
-  a 

by  (5.2a,  b,  c)  with  k  >  r. 

Then  there  exist  an  hQ  independent  of  e,  such  that  for  all  h  <  hQ,  there 

exists  a  unique  solution  u  £  S  to 

u  r 


(k) 


(5.9a)  B  (u  ,v  )  =  F (v  )  for  each  v  £S 

L  a  a  a  a  a 


Also,  let  u  £S  be  the  unique  solution  to 


(5.9b)  B^(uL,vL)  *  F (v^)  for  each 
Then,  for  1  ^  p  <  00  , 

JL  ”  x 1  1 H 


p,A 


r  ,k+ ^  1  I  r  (k)  I  ] 

C1  max  h  ||fo  !Il(i) 
2  00  J 


+  C2  max  hk+1  r{  |  | a (k+1_r) | |  r .  .  +  j  |b 

j  r 


l  lu(k-r) 


Hr  (I  . ) 

00  .1 


} 


with  Cj  and  C2  independent  of  e  and  A. 


4  b 


J 


Proof .  The  existence  and  uniqueness  of  u  will  be  established  if  the  homogeneous 


problem, 

(5.10) 


B.(w,v  )  =  0  for  each  v  C  S^ 

A  a  a  a 


has  onlv  the  zero  solution  in  S  . 

Suppose  (5.10)  holds.  Let  v.  C  be  s.t.  v.(x.)  =  v  (x.),  j  =  0,.. 

L  L  J  *  1 

Then  =  P^(v^).  From  (5.10)  it  follows  that 

BA(w’vL)  =  BA(w’VVa}  for  each  \C  SLn)’  and  vct  C  ^  ’ 

and  hence 


(5.11) 


ba(w*V 


!ba(w-VV 


l.e.A 


l,e,A 


for  each  vL  £  and  vn  £S^. 


Because  the  inf-sup  condition  holds  for  B^(-,.)  over  the  spaces  Sf  x  S^  for  h 
sufficiently  small  (Theorem  3.2  &  Theorem  4.3),  there  exists  a  such 

that  the  left  hand  side  of  (5.11)  is  larger  than  c||wj|  o  .  The  right  hand 

°°,  A 

side  of  (5.11)  can  be  bounded  by  Lemma  5.2,  and  hence  it  follows  that 


(5.12) 


CjJ  I wf  |Ho  <  C2hk+1"r| |w| |Ho 

°°,A  °°,A 


where  and  C2  are  independent  of  h  and  £.  Since  k  r,  there  exists  an  h^  such 

that  for  h  £  h  ,  the  only  way  that  (5.12)  can  be  satisfied  is  if  | |w| |^o  =  0, 

—  0  ®,A 


which  implies  w  =  0. 


Since 


it  follows  from  (5.9a,b)  that 


YWV  =  F(VV«)  -  WW' 


From  Theorem  3.2  it  follows  that 


(5.13) 


F  (vL-vf,), 


'UL  ”  Ua  '  o  1  C1  SUP  -1 - — 

Hp,A  vCS(n)  MVlII  2 

L  L  ri 


q  >e  .A 


+  sup 


vLes^>  ||vl||h2 


for  each  v  and  1  <  p  <  where  —  +  —  =  1. 

ec  tt  ’  —  1  —  p  q 

I  F  (v  -v  ) I 

First,  consider  the  term  sup  '  L  ,x  l 

vCS<n)  1  Ti-I  |„2 

L^bL  p,e,A 


q  ,  C  ,  A 


Let  F . (x) 
J 


-A, 


f  (s)dx 
o 


for  each  x  C  I..  Let  F^  be  the  polynomial  of  degree  k  on  1^  such  that 

F,(x.+ih./k)  =  F.(x.+ih./k)  for  i  =  0,1,  ...,k.  Then 
.1  J  .1  J  J  1 


(5.14) 


'FJ  "  FiML  (I.)  < 


I  lr  (k)  I  I  ,  k+1 

llfo 

_ 00  J 

4 (k+1) ! 


As  before,  let  v  £  S'  be  such  that  v  (x.)  =  vT(x.),  i  =  0,...,N.  Then, 

a  a  a  j  L  3 


F  (VV 


N  ^ 

I  f  (x) (v  -V  ) (x)dx 
.1  =  1  */1 .  0  n 


N  r 

-  Z  /f, 

1=1  I.  1 
-1 


(x  )  (Vj -v  )  (x)dx 
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N 

=  L  /  (F  -F  )(v'-V)dx, 
.1  =  1  1  1  L  a 


the  last  equality  following  from  (5.3).  Thus, 


(5.15) 


1f(v  -V  )  i  <  t  ! IF1  -  ^llL  (i.)Mv’  -  v;l!  (  , 

L,  J  =  1  °°  V  1  j 


From  (5.15),  (5.14),  Lemma  5.1,  inequality  (2.18),  and  the  fact  that 


I  v  j  |  2  <  |  jv'  |  2  ,  1  <  q  <  °°,  it  follow's  that 

H  7  A  —  H  .  —  — 

1 ;  e,A  q  ,e,A 


(5.16) 


VLCSf 


sup 

-<,(n) 


F  (vL-V()  | 


r«-N„j 


—  C  max  h™  (I  V 

.1  "  .1 


q.e.A 


Using  Lemma  5.2,  and  inequality  (5.16)  it  follows  from  (5.13)  that 


(5.17) 


!uj-u  I  |  <  C1  max  h^4^  |  |f  ^  I  |  .  . 

“L  cr'Ho  -  1  .1  11  o  'L JJ  ) 

p,  A 


,  r  ,  k+l-r r |  |  (k+l-r) i  i 

+  C2  max  h.  { |  |a  |  lHr(].  ) 


+  lh<k'r)U„r,T  s> 


Hr(I.) ' 1 |Uu‘ 'L  (I.)  . 

oo  .1  oo  J 


Bv  hypothesis,  |  |u|  is  bounded  indenendentlv  of  G.  Since  Uj  is  a  quasi- 

00 

optimal  anproximation  to  u,  it  follows  that  I  | uT  |  !  t  is  bounded  independently  of 

"oo 

G.  From  (5.17)  it  follows  that  if  h  is  sifficiently  small,  then  | |u  | j  is  also 

^  Loo 

bounded  independently  of  e.  That  J[u  ||.  is  bounded  independently  of  G,  and  A 

a  Loo 

(provided  h  is  sufficiently  small),  combined  with  (5.17),  proves  the  theorem. 

This  theorem  shows  that  if  we  project  the  space  onto  k  upwinded 

(k)  (k) 

polynomials  in  ,  then  the  finite  element  solution  witn  as  the  test  snace 

will  have  an  error  composed  of  two  parts: 
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1  -  <  1  !u  -  u,  I !  +  ;  |u  -  u  1 1 


'!  •  “*  A~  ”l.-o 


P.A 


1  '*  H° 


T*  .  first  part 


is  quasi -opt  ima!  and  hence  the  'nest  order  of  thi- 


p ,  A 


r+1 

term  that  can  he  expected  is  0(h  ).  Tf  k  =  2r  +  1,  the  theorem  savs  that  the 


second  term  Mu,  -  u  M 

1.  v  Ho 


0(h  '  ),  which,  is  one  higher  order  than  the 


P.A 


optimal  error. 


Coroll a rv  5.4.  Suppose  that  all  the  assumnt ions  of  Theorem  5.3  hold  with  k=2r+l . 
Then  there  exists  an  hD  independent  of  r  such  that  for  all  h  <  hQ 

I iu  "  u  ! IH°  5  Ci  lnf  llu  -  w| lHo  +  C2  hr+“  . 
p,A  p.A 

r+° 

We  remark  that  when  using  a  polynomial] v  upwinded  test  space  the  term  C0h 

is  unadvoidable .  It  was  proven  in  [3]  and  [34]  that  quasi-opt imal i t v  is 

unobtainable  when  a  test  space  containing  basis  functions  of  the  form  '5.1),  with 

g^(x)  independent  of  e  and  h,  is  used.  However,  if  local  smoothness  on  a,  •>  mil 

f  is  assumed-a  condition  alwavs  satisfied  is  practice,  we  can  obtain  an  additional 
y-f-2 

error  of  0(h  ).  This  error  can  in  general  be  neglected  because  tin'-  best 

r+  ] 

approximation  from  a  function  w£7>  has  order  0(h  ). 

Because  of  the  second  part  of  this  work  -  the  a-posteriori  error  estimates, 

it  is  important  to  keep  the  projection  error  of  one  higher  order  than  the  optimal 

error.  This  should  also  be  true  of  the  numerical  quadrature  errors.  A 

quadrature  rule,  which  is  exact  when  n(x)  is  a  piecewise  polynomial  of  degree  r+2  , 

and  h  and  f  are  piecewise  polynomials  of  degree  r+1,  is  derived  in  [34],  and  shown 

t+2 

to  produce  an  error  of  order  0(h  )  as  well. 


->0 


UNCLASSIFIED 


CHAPTER  A 
CONCLUSION 

In  this  paper  it  has  been  shown  that  quisi-opt imality  is  obtainable  for  a 
finite  element  solution  when  the  test  space  is  composed  of  functions  which  are 
"nearly"  L*  splines.  The  norm  used  to  measure  the  errors  is  very  close  to  an 
Lp  norm  which  is  important  particularlv  if  the  location  and  shape  of  the  boundary 
layer  are  important.  Furthermore,  although  it  was  shown  in  [3]  and  [34]  that 
quasi-optimalitv  is  unobtainable  when  using  a  polvnomiallv  upwinded  test  space, 
we  have  shown  that  a  "nearly"  quasi-opt imal  result  is  possible  if  the  input 
functions  a(x),  b(x)  and  f (x)  are  piecewise  smooth.  This  "nearly"  quasi-opt imal 
result  is  sufficient  for  finding  a-posteriori  error  estimates  and  proving  that  the 
error  estimate  converges  to  the  true  error  as  h  =  h(A)  -  0.  This  result  is  proven 
and  an  adaptive  mesh  refinement  procedure  and  numerical  results  are  presented 
in  the  second  part  of  this  paper  [35]. 

Many  of  the  results  of  this  paper  (particularly  the  embedding  result)  used 
bounds  on  the  Green's  function.  For  turning  point  problems,  in  which  there  is  a 
point  xG  C  I  such  that  a(xQ)  =  0,  the  Green's  function  is  not  bounded  independently 
of  C.  However,  in  [  6  ]  sharp  bounds  are  given  for  the  Green's  functions  of 
turning  point  problems.  The  numerical  results  for  a  turning  point  problem  are 
given  in  Part  II,  [35]  and  suggest  that  analogous  results  hold  for  this  case 
as  well. 

Upwinding  can  easily  be  implemented  in  two  dimensional  problems.  On  a 
rectangular  mesh,  upwinding  can  be  done  bv  simply  upwinding  in  each  direction 
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separately.  This  is  done  in  [19]  -  [21]  for  >-  quadratic  unwinding,  in  [23]  for 

upwinding  the  integration  point,  and  in  [21]  for  tipwinded  elements  which  are  local 

asymptotic  expansions  of  the  solution. 

In  [7],  Brooks  and  Hughes  describe  a  method  in  which  upwinding  is  done  only 

in  the  direction  of  the  flow  ("streamline  diffusion  method).  A  mathematical  analv- 

sis  of  this  method  was  performed  by  Navert  [31].  Interior  error  estimates  were 

k+h 

proven  to  decrease  with  rate  h  in  which  is  h  of  a  rower  lower  than  the 
optimal  rate.  However,  the  optimal  rate  was  observed  in  the  numerical  results. 

The  polynomial  upwinding  presented  in  this  p^per  can  also  be  implemented  in 
two  dimensions.  Nevertheless,  major  theoretical  questions  still  remain  unsolved. 


32 


REFERENCES 


1.  I.  Babuska,  A.  K.  Aziz,  Survey  Lectures  on  the  Mathematical  Foundation 

of  the  Finite  Element  Method,  A.  K.  Aziz,  ed. ,  The  Mathematical  Foundations 
of  the  Finite  Element  Method  with  Applications  to  Partial  Differential 
Equations,  Academic  Press,  New  York,  1973. 

2.  I.  Babuska,  J.  Osborn,  Analysis  of  Finite  Element  Methods  for  Second 
Or  ier  Boundary  Value  Problems  Using  Mesh  Dependent  Norms,  Numerische 
Mathematik,  Vol.  34,  1980,  pp.  41-62. 

3.  I.  Babuska,  W.  G.  Szvmczak,  An  Error  Analysis  for  the  Finite  Element 
Method  Applied  to  Convection  Diffusion  Problems,  Comp.  Meth.  Appl,  Mech. 
and  Engng. ,  Vol.  31,  (1982),  pp.  19-42. 

4.  J.  W. ,  Barrett,  K.  W.  Morton,  Optimal  Petrov-Galerkin  Methods  Through 
Approximate  Symmetrization,  IMA  Journal  of  Numerical  Analysis,  Vol.  I, 
(1981),  pp.  439-468. 

5.  _ ,  Optimal  Finite  Element  Solutions  to  Diffusion  Convection 

Problems  in  One  Dimension,  Int.  J.  Num.  Meth.  Engng.,  Vol.  15,  (1980). 

6.  A.  E.  Berger,  H.  Han,  R.  B.  Kellogg,  A  Priori  Estimates  and  Analysis 
of  a  Numerical  Method  for  a  Turning  Point  Problem,  to  appear. 

7.  A.  Brooks,  T.  J.  R.  Hughes,  Stream! ine-Upwind/Petrov-Galer kin  Methods 
for  Advection  Dominated  Flows,  Third  Internationa]  Conference  on  Finite 


Element  Methods  in  Fluid  Flow,  1980. 


r 


1 


j 

i 

k 

1 

"i 

j»« 

* 

I 


8.  J.  Christie,  D.  F.  Griffiths,  A.  R.  Mitchell,  0.  C.  Zienkiewicz, 

Finite  Element  Methods  for  Second  Order  Differential  Equations  with 
Significant  First  Derivatives,  Int .  J.  Num.  Meths.  in  Engng.,  Vol.  10, 
(1976),  pp.  1389-1396. 

9.  J.  Christie,  A.  R. ,  Mitchell,  Upwinding  of  High  Order  Galerkin  Methods 

in  Conduction-Convection  Problems,  Int.  J.  Num.  Meths.  in  Engng.,  Vol.  32, 
(1978),  pp.  1764-1771. 

10.  J.  D.  Cole,  J.  Kevorkian,  Perturbation  Methods  and  Applied  Mathematics, 
Springer-Verlag,  New  York,  1981. 

11.  P.  P.  N.  De  Groen,  A  Finite  Element  Method  With  Large  Mesh  Width  For 
Stiff  Two-point  Boundary  Value  Problems,  Preprint,  Eindhoven,  Univ. , 
Department  of  Mathematics,  Eindhoven,  the  Netherlands,  1978. 

12.  P.  P.  N.  De  Groen,  P.  W.  Hemker,  Error  Bounds  for  Exponentially  Fitted 
Gerlerkin  Methods  Applied  to  Stiff  Two-point  Boundary  Value  Problems, 

P.  W.  Hemker,  J.  J.  H.  Miller,  ed..  Numerical  Analysis  of  Singular 
Perturbation  Problems,  Academic  Press,  New  York,  1979. 

13.  R.  F.  Diaz-Munio,  C.  Well  ford,  Jr.,  A  Finite  Element  Singlular  Perturbation 
Technique  for  Convection-Diffusion  Problems.  Part  I:  The  One-Dimensional 
Case,  A.S.M.E.  Journal  of  Applied  Mechanics,  to  appear. 

14.  _ ,  A  Finite  Element-Singular  Perturbation  Technique  for 

Convection-Diffusion  Problems,  Part  II:  Two  Dimensional  Problems, 

A.S.M.E.  Journal  of  Applied  Mechanics,  to  appear. 


J 


15.  W.  Eckhaus,  Asymptotic  Analysis  of  Singular  Perturbations,  North-Hollar'1 
Pub.  Co.,  Amsterdam,  New  York,  1979. 

16.  P.  M.  Gresho,  R.  L.  Lee,  Don't  Suppress  the  Wiggles-They ' re  Telling  You 
Something!,  Computers  and  Fluids,  Vol.  9,  (1981),  pp.  223-253. 

17.  D.  F.  Griffiths,  A.  F.  Mitchell,  On  Generating  Upwind  Finite  Element 
Methods,  T.  J.  Hughes,  ed.,  Finite  Elements  for  Convection  Dominated 
Flows,  (AMD,  1979),  pp.  91-104. 

18.  D.  F.  Griffiths,  J.  Lorenz,  An  Analysis  of  the  Petrov-Galerkin  Finite 
Element  Method,  Comp.Meth.  Appl.  Mech.  Engng. ,  Vol.  14,  (1978),  pp.  39-64. 

19.  J.  C.,  Heinrich,  P.  S.  Huyakom,  0.  C.  Zienkiewicz,  A.  R.  Mitchell, 

An  'Upwind'  Finite  Element  Scheme  for  Two-Dimensional  Convective  Transport 
Equation,  Int .  J.  Num.  Meth.  Engng.,  Vol.  11,  (1977),  pp.  131-143. 

20.  J.  C.  Heinrich,  0.  C.  Zienkiewicz,  Quadratic  Finite  Element  Schemes  for 
Two-Dimensional  Convective-transport  Problems,  Int.  J.  Num.  Meth.  Engng., 

Vol.  11,  (1977),  1831-1844. 

21.  _ ,  The  Finite  Element  Method  and  Upwinding  Techniques  in  the 

Numerical  Solution  of  Convection  Dominated  Flow  Problems,  T.  J.  R.  Hughes,  ed. 
Finite  Element  Methods  for  Convection  Dominated  Flows  (AMD,  1979),  pp. 

105-135. 

22.  P.  W.  Hemker,  A  Numerical  Study  of  Stiff  Two-Point  Boundary  Problems, 
Mathematisch  Centrum,  Amsterdam,  1977. 

23.  T.  J.  R.  Hughes,  W.  K.  Liu,  A.  Brooks,  Finite  Element  Analysis  of 
Incompressible  Viscous  Flows  by  the  Penalty  Function  Formulation,  J.  Comp. 
Physics,  Vol.  30,  (1979),  pp.  1-60. 


C.  Johnson,  U.  Navert,  An  Analysis  of  Some  Finite  Element  Methods  for 
Convection-Diffusion  Problems,  Analytical  and  Numerical  Approaches  to 
Asymptotic  Problems  in  Analysis,  North  Holland,  1981. 

C.  Johnson,  Finite  Element  Methods  For  Convection  Diffusion  Problems, 
Technical  Report,  Chalmers  Univ.  of  Technology,  Goteberg,  Sweden  1981. 

R.  B.  Kellogg,  H.  Han,  Numerical  Analysis  of  Singular  Perturbation 
Problems,  Univ.  of  Maryland,  I.P.S.T.,  Technical  Note  BN-977,  1981. 

R.  B.  Kellogg,  H.  Han,  The  Finite  Element  Method  For  A  Singular  Perturbation 
Problem  Using  Enriched  Subspaces,  Univ.  of  Maryland,  I.P.S.T.,  Technical 
Note  BN-978,  1981. 

R.  B.  Kellogg,  A.  Tsan,  Analysis  of  Some  Difference  Approximations  for 
A  Singular  Perturbation  Problem  Without  Turning  Points,  Math,  of  Comp., 

Vol .  32,  (1978),  pp.  1025-1039. 

•  • 

J.  Lorenz,  Zur  Theorie  und  Numerik  von  Dif ferenzenverfahren  fur  Singulare 

M  •• 

Storungen,  Universitat  Konstanz,  1981. 

A.  Mitchell,  D.  F.  Griffiths,  Upwinding  by  Petrov-Galerkin  Methods  in 
Convection-Diffusion  Problems,  University  of  Dundee  Department  of  Mathe¬ 
matics  Report  NA/36,  1979. 

U.  Navert,  A  Finite  Element  Method  for  Convection-Diffusion  Problems, 

Thesis, Chalmers  Univ.  of  Technology,  1982. 

M.  H.  Protter,  H.  J.  Weinberger,  Maximum  Principles  in  Differential 
Equations,  Prentice  Hall,  Inc.,  Englewood  Cliffs,  New  Jersey,  1967. 


33.  I.  Stakgold,  Boundary  Value  Problems  of  Mathematical  Physics,  Vol.  1, 

The  Macmillan  Company,  New  York,  1967. 

34.  W.  G.  Szymczak,  An  Adaptive  Finite  Element  Method  for  Convection-Diffusion 
Problems,  Thesis,  Univ.  of  Maryland,  1982. 

35.  W.  G.  Szymczak,  I.  Babuska,  An  Analysis  of  A  Finite  Element  Method  For 
Convection-Diffusion  Problems.  Part  II:  A-Posteriori  Error  Estimates 
and  Adaptivity,  to  appear. 

36.  M.  I.  Visik,  L.  A.  Lyusternik,  Regular  Degeneration  and  Boundary  Layer 
For  Linear  Differential  Equations  With  Small  Parameter,  Amer.  Math.  Soc., 
trans.  Vol.  20,  No.  2,  (1962),  pp.  239-364. 

37.  0.  C.  Zienkiewicz,  J.  C.  Heinrich,  The  Finite  Element  Method  and 
Convection  Problems  in  Fluid  Mechanics,  Gallagher,  Zienkiewicz,  Oden 
Morandi,  Cecchi,  and  Taylor,  eds. ,  Finite  Elements  in  Fluids,  Wiley, 
London,  1978,  pp.  1-21. 


